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Conversion  Factors, 

Non-Si  to  SI  Units  of  Measurement 

Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units  as  follows: 


Multiply 

By 

To  Obtain 

acres 

4,046.873 

square  meters 

degrees  (angle) 

0.01745329 

radian 

feet 

0.3048 

meters 

knots  (international) 

0.5144444 

meters  per  second 

miles  (US  Statute) 

k6093 

kilometers 

nautical  miles 


1.852 


kilometers 


1  Introduction 


Background 

Kikiaola  Harbor  is  a  small,  shallow-draft  harbor  on  the  Island  of  Kauai.  The  harbor 
is  approximately  100  miles  west  and  a  little  north  of  Honolulu,  Oahu,  and  is  located 
along  the  western  part  of  Kauai’s  south  shore  (Figure  1).  Nearby  towns  of  Kekaha  and 
Waimea  are  1  mile  northwest  and  1.5  mile  east  of  the  harbor,  respectively.  Lihue,  the 
county  seat  and  business  center  of  Kauai,  is  located  approximately  23  miles  east  of  the 
harbor.  The  local  shoreline  is  a  relatively  straight,  low,  wide  beach  that  reaches  from 
Oomano  Point  at  its  western  extent  to  the  Waimea  River  on  the  east,  a  distance  of  2.7 
miles. 

Kikiaola  Harbor  was  originally  developed  by  die  State  of  Hawaii  in  1959. 

Inner  and  outer  stub  extensions  to  the  east  breakwater  and  a  short  inner 
breakwater  were  added  in  1964  to  form  the  present  harbor  (Figure  2).  The 
additional  structures  were  needed  to  reduce  surge  within  the  harbor.  A  wharf  and 
boat  ramp  are  located  along  die  north  boundary  of  the  harbor,  east  of  the  inner 
breakwater. 

Prevailing  northeast  tradewinds  result  in  a  strong  predominance  of  winds  from 
northeast,  east,  and  southeast  at  Kikiaola  Harbor.  Typical  wind  speeds  are  10  to 
20  mph.  Winter  storms  can  generate  strong  winds  from  die  south.  The  harbor  is 
exposed  to  waves  approaching  from  a  sector  between  N  82°  W  and  S  46°  E 
(Figure  1),  though  the  small  island  of  Niihau  creates  some  sheltering  in  die 
western  part  of  this  sector.  Southern  swell,  generated  by  storms  in  die  southern 
Pacific  and  Indian  Oceans,  is  a  significant  part  of  the  wave  climate.  Also,  waves 
generated  by  storms  in  the  North  Pacific  can  wrap  around  the  western  side  of 
Kauai  and  affect  Kikiaola  Harbor.  Hurricanes  can  attack  the  harbor.  This  source 
of  waves  is  important  for  structure  design  but  is  sufficiendy  rare  that  it  does  not 
impact  the  operational  concerns  of  die  present  study. 

Use  of  the  existing  harbor  is  limited  by  two  primary  factors.  First,  the  harbor 
is  quite  shallow.  Sediment  movement  along  die  local  coast,  predominandy  from 
east  to  west,  has  resulted  in  shoaling  of  die  entrance  and  inner  harbor.  Second, 
die  existing  entrance  experiences  breaking  wave  conditions  which  are  hazardous 
to  navigation.  These  two  factors  are  interrelated.  Breaking  waves  are  more 
likely  in  die  existing,  shoaled  entrance  than  they  would  be  in  a  deeper,  maintain^ 
entrance  channel. 
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Figure  1.  Kikiaola  Harbor  location 


The  U.S.  Army  Engineer  Division,  Pacific  Ocean,  developed  several  plans  for 
improving  operational  conditions  in  Kikiaola  Harbor  (U.S.  Army  Engineer 
Division,  Pacific  Ocean  1980).  The  present  study  focuses  on  evaluating  wave 
conditions  in  the  most  promising  plans  to  help  insure  that  operational  objectives 
are  met. 

Study  Approach 

The  study  described  in  this  report  was  performed  by  die  U.S.  Army  Engineer 
Waterways  Experiment  Station  (WES),  Coastal  and  Hydraulics  Laboratory 
(CHL).  The  approach  consisted  of  the  following  components: 

a.  Develop  wave  climate  information  at  the  harbor  site. 

b.  Use  a  numerical  model  to  investigate  existing  and  alternative  harbor 
modification  plans. 
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Preface 


This  study  was  authorized  by  the  U.S.  Army  Engineer  Division,  Pacific  Ocean 
(POD),  and  was  conducted  by  personnel  of  the  Coastal  Hydrodynamics  Branch 
(CHB),  Navigation  and  Harbors  Division  (NHD),  Coastal  and  Hydraulics 
Laboratory  (CHL),  of  the  U.S.  Army  Engineer  Waterways  Experiment  Station 
(WES).  The  study  was  conducted  during  the  period  August  1996  through  July 
1997.  Messrs.  Tim  Young  and  Lincoln  Gayagas,  POD,  oversaw  progress  of  the 
study. 

Dr.  Edward  F.  Thompson,  CHB,  was  the  WES  point  of  contact  for  the  study. 
This  report  was  prepared  by  Dr.  Thompson,  Dr.  Lihwa  Lin,  and  Ms.  Lori  L. 
Hadley,  all  of  CHB,  and  Dr.  Jon  M.  Hubertz,  formerly  of  CHB.  Drs.  Lin  and 
Hubertz  conducted  the  wave  climate  analysis  portion  of  the  study.  Direct  supervi¬ 
sion  was  provided  by  Dr.  Martin  C.  Miller,  former  Chief,  CHB.  General 
supervision  was  provided  by  Mr.  C.  E.  Chatham,  Chief,  NHD,  Mr.  Charles  C. 
Calhoun,  Jr.,  Assistant  Director,  CHL,  and  Dr.  James  R.  Houston,  Director,  CHL. 

At  the  time  of  publication  of  this  report.  Dr.  Robert  W.  Whalin  was  Director  of 
WES.  COL  Robin  R.  Cababa,  EN,  was  Commander. 
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Figure  2.  Kikiaofa  Harbor,  existing  plan 


Wind  wave  and  swell  climate  was  investigated  primarily  with  numerical 
hindcast  information  covering  a  period  of  one  year.  Ideally,  a  longer  time  period 
would  be  used  to  establish  wave  climate,  but  only  one  year  was  available  from  a 
comprehensive  re-hindcast  of  the  Pacific  Ocean  within  the  time  frame  of  this 
study.  Even  one  year  of  information  can  be  expected  to  give  reasonable 
representation  of  die  lower  99  percent  of  wave  conditions,  which  cover  die  main 
concerns  in  this  study  (Thompson  and  Harris  1972).  Buoy  measurements  from 
several  locations  were  used  to  help  validate  the  hindcasts.  Hindcast  wave 
information  was  used  as  a  boundary  condition  for  nested  finer  grids  which  allowed 
sheltering  of  the  islands  of  Kauai  and  Niihau  to  be  modeled.  Deepwater  waves 
offshore  of  Kikiaola  Harbor  were  transformed  by  a  simple  method  to  a  depth  of  4 
m  for  use  as  an  incident  condition  in  die  harbor  wave  model.  The  wave  climate 
study  is  presented  in  Chapter  2. 

A  numerical  wave  model  was  set  up  to  cover  die  harbor  and  an  area  outside 
die  harbor  extending  about  300  yd  seaward  of  die  entrance.  Two  proposed  harbor 
plans  and  the  existing  harbor  were  studied.  Both  plans  include  modification  of  die 
east  and  west  breakwaters,  dredging  of  an  entrance  channel  to  a  depth  of  12  ft 
below  MLLW,  and  dredging  of  an  inner  access  channel  to  a  depth  of  8  ft 
MLLW.  Special  features  of  each  plan  are: 

a.  Plan  1  (Figure  3).  Remove  outer  stub  of  east  breakwater;  remove  and 
reconstruct  inner  stub  of  east  breakwater  a  small  distance  further  east;  raise 
crest  elevation  of  exposed  portions  of  east  breakwater  by  3-4  ft  and  flatten 
seaward  slope  to  1:2;  widen  outer  220  ft  of  west  breakwater;  dredge  725-ft 
long  entrance  channel  with  width  varying  from  10S  to  205  ft  and  maneuvering 
area  to  facilitate  a  90  deg  turn  into  access  channel;  dredge  320-ft  long  access 
channel  varying  in  width  from  70  to  105  ft. 

b.  Plan  6  (Figure  4).  Remove  outer  and  inner  stubs  of  east  breakwater;  raise 
crest  elevation  of  exposed  portions  of  east  breakwater  by  3-4  ft  and  flatten 
seaward  slope  to  1:2;  extend  east  breakwater  further  west  to  a  distance  of  100 
ft  past  the  existing  west  breakwater  location;  shorten  west  breakwater  to  allow 
space  for  access  channel;  dredge  entrance  and  access  channels  comparable  to 
those  in  Plan  1. 

The  numerical  model  used  for  the  studies,  HARBD,  is  the  standard 
USAEWES  tool  for  numerical  harbor  wave  investigations.  The  model  includes 
the  following  assumptions: 

a.  No  wave  transmission  through  the  breakwaters. 

b.  No  wave  overtopping  of  structures. 

c.  Structure  crest  elevations  above  the  water  surface  cannot  be  tested  or 
optimized. 

d.  Currents  in  tire  channel  can  not  be  evaluated. 
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Figure  4.  Plan  6 


e.  Wave  breaking  effects  in  the  entrance  and  harbor  cannot  be  considered 
directly. 

/  No  nonlinear  effects  are  considered. 

g.  Diffraction  around  structure  ends  is  represented  by  diffraction  around  a 
blunt  vertical  wall  with  specified  reflection  coefficient. 

Despite  limitations  imposed  by  die  above  assumptions,  HARBD  is  considered 
suitable  for  meeting  die  numerical  modeling  objectives  of  the  Kikiaola  Harbor 
study.  Development  of  the  numerical  model  and  test  procedures  is  described  in 
Chapter  3. 

Response  of  the  existing  and  alternative  harbor  plans  to  waves  was  studied 
using  numerical  model  results.  Harbor  response  to  wind  waves  and  swell  (short 
waves )  is  presented  in  Chapter  4.  The  harbor  short  wave  response  is  related  to 
wave  climate  and  to  standard  U.S.  Army  Corps  of  Engineers  criteria  in  channels 
and  berthing  areas.  Harbor  oscillation  characteristics  (response  to  long  waves) 
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are  presented  in  Chapter  5.  The  long  wave  study  included  only  die  existing 
harbor  and  Plan  6. 

Conclusions  and  recommendations  are  given  in  Chapter  6.  This  chapter  is 
followed  by  references  and  appendices  with  detailed  information  supporting  the 
main  report  and  notation  definitions. 
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2  Wind  Wave  and  Swell 
Climate 


Sources 


Four  sources  of  wind  wave  and  swell  information  were  available  to  develop 
wave  climate  outside  die  harbor  entrance,  including  three  National  Data  Buoy 
Center  (NDBC)  buoys  with  open  exposure  to  die  south  and  the  Wave  Information 
Studies  (WIS).  WIS  has  hindcast  waves  over  die  Pacific  Ocean  and  saved 
information  at  selected  deepwater  stations  around  die  Hawaiian  Islands.  Buoy 
locations  and  corresponding  WIS  deepwater  (Level  1)  stations  are  shown  in  Fig. 

5.  The  original  WIS  Pacific  hindcast  covered  only  die  north  Pacific  (Corson  et 
al.  1986).  That  study  is  presendy  being  updated  to  extend  coverage  into  die  south 
Pacific  as  well.  At  die  time  of  the  Kikiaola  Harbor  study,  only  one  year  (1989) 
was  completed.  Because  waves  from  die  south  Pacific  are  a  critical  part  of  die 
climate  at  Kikiaola  Harbor,  the  1-year  updated  WIS  hindcast  was  used  in  this 
study  in  preference  to  die  original  WIS  information. 

Deepwater  Wave  Climate 

The  deepwater  WIS  hindcast  for  1989  was  calculated  in  three  steps,  with  each 
step  giving  increased  refinement.  The  initial  and  coarsest  step,  Level  1,  covered 
the  entire  Pacific  Ocean  basin  with  a  grid  spacing  of  2.5  deg  of  latitude/longitude 
(Fig.  6).  Information  at  a  sequence  of  points  enclosing  U.S.  Pacific  coasts 
(shown  in  the  figure)  was  saved  to  use  as  a  boundary  condition  for  a  more 
refined,  localized  Level  2  grid.  The  Level  2  grid  (Fig.  7),  with  mesh  spacing  1.0 
deg  latitude/longitude,  provided  boundary  conditions  around  the  Hawaiian  Islands 
for  a  Level  3  grid  with  0.25-deg  resolution  (Fig.  8).  An  additional  refinement 
was  added  for  this  study  to  represent  the  sheltering  effects  of  the  islands  of  Kauai 
and  Niihau  (the  small  island  15  miles  southwest  of  Kauai)  on  the  project  site. 

This  Level  4  grid  had  a  resolution  of  0.017-deg  latitude/longitude  and 
encompassed  both  Kauai  and  Niihau.  Wind  wave  growth  was  included  in  Levels 
1-3  and  propagation  effects  were  included  in  all  levels  of  hindcast. 
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Figure  5.  Location  map  for  NDBC  buoys  and  nearby  WIS  Level  1  stations 


Hindcast  waves  during  1989  were  compared  with  NDBC  buoy  measurements 
over  the  same  time  period.  Each  buoy  was  compared  to  the  nearest  WIS  Level  1 
smtion.  Comparison  statistics  show  root-mean-square  (rms)  differences  and  mean 
differences  in  significant  wave  height,  Hs,  and  mean  spectral  wave  period,  Tm 
(Table  1).  The  biases  are  quite  small  and  die  rms  differences  are  typical  of  a 
validated  hindcast  model  (Brooks  and  Brandon  1995). 

It  is  important  to  evaluate  whether  die  time  period  of  special  hindcast  is 
representative  of  the  long  term  climate  incident  to  the  south  coast  of  Kauai.  Wave 
parameter  summaries  for  the  deepwater  sources  are  compared  in  Table  2.  Long 
term  mean  values  of  H,  for  the  three  buoys  are  within  0. 1  m  of  the  corresponding 
WIS  stations  during  1989.  Mean  Tp  values  are  within  1  sec.  Standard  deviations 
of  Hs  are  similar  between  the  buoys  and  WIS  stations.  Standard  deviation  of  Tp  is 
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Figure  6.  Location  map  for  wave  climate  study,  WIS  Level  1  grid 


Table  1 

Comparison  of  WIS  Level  1  Hindcasts  and  NDBC  Buoy  Data,  1989 


WIS 

Station 

NDBC 

Buoy 

RMS  Difference 

Bias  I 

H.{ m) 

r.(sec) 

HA m) 

7.  (sec)  | 

L1.1 

51001 

0.47 

1.6 

-006 

-0.6  \ 

LI  .2 

51002 

0.53 

1.9 

0.03 

-04  | 

LI  .3 

51003 

0.51 

1.6 

0.10 

-0.6  || 

'  Calculation  based  on  model  minus  measured  values 
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Figure  7.  Location  map  for  wave  climate  study,  WIS  Level  2  grid 


somewhat  greater  for  the  buoys  than  for  die  WIS  stations.  Overall,  die  WIS 
information  for  die  year  1989  appears  representative  of  die  long  term  climate. 


Wave  information  from  die  special  WIS  Level  4  hindcast  was  saved  at  a 
deepwater  point  about  one  mile  offshore  from  Kikiaola  Harbor.  Wave  climate  at 
this  point  is  summarized  in  Figs.  9-11.  Wave  directions  are  predominandy  from 
the  south  (180  deg)  and  west  northwest  (300  deg).  Waves  coming  from  the  west, 
northwest,  and  north  outside  die  islands  are  partially  blocked  by  Niihau  and 
western  Kauai  before  they  can  arrive  at  Kikiaola  Harbor.  The  concentration  of 
waves  coming  from  300  deg  indicates  that  a  significant  amount  of  wave  energy  is 
penetrating  around  and  between  die  islands. 
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Figure  8.  Location  map  for  wave  climate  study,  WIS  Level  3  grid 


Table  2 

Summary  Statistics,  NDBC  Buoys  and  WIS 

Statistical  Parameter 

Buoy 

610011 

WIS 

LI. I3 

Buoy 

61002* 

WIS 

L12' 

Buoy 

61003* 

WIS 

L1.3* 

Mean  H,  (m) 

2.4  * 

2.3 

24 

2.3 

2.2 

2.3 

Standard  deviation  of  H,  (m) 

0.9 

0.6 

07 

0.5 

0.7 

0.6 

Mean  7,  (sec) 

105 

11 

10.2 

11 

10.4 

11 

Standard  deviation  of  7,  (sec) 

2.9 

2 

2.8 

2 

2.9 

1 

1  Data  from  Feb  81  through  Dec  93 

3  Data  from  Jan- Dec  89 

3  Data  from  Sep  84  through  Dec  93 

4  Data  from  Nov  84  through  Dec  93 
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Figure  9.  Deepwater  wave  climate,  Ht 


Figure  1 0.  Deepwater  wave  climate,  Tp 


Figure  1 1 .  Deepwater  wave  climate,  6m  (deg,  coming  from) 
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Wave  Climate  at  Kikiaola  Harbor 


The  deepwater  wave  climate  analysis  suggests  that  data  from  die  special  one- 
year  hindcast  reasonably  characterizes  die  wave  climate  immediately  incident  to 
Kikiaola  Harbor.  The  hindcast  information  must  be  transformed  into  shallow 
water  to  provide  wave  climate  at  the  seaward  boundary  of  die  HARBD  model,  in 
water  depth  of  approximately  13  ft  (4  m). 

Initially,  die  nearshore  transformation  was  attempted  using  die  STWAVE  and 
RCPWAVE  models  in  die  Automated  Coastal  Engineering  System  (ACES  2.0) 
(Leenknecht  and  Tanner  1996;  Leenknecht,  Tanner,  and  Sherlock  1997).  The 
models  gave  inconsistent  results  over  die  highly  irregular  bathymetry  near 
Kikiaola  Harbor.  Since  bathymetry  contours  are  fairly  straight  and  parallel 
seaward  of  about  18  ft  MLLW,  die  deepwater  wave  conditions  were  transformed 
to  13-ft  depth  using  simplified  refraction  procedures.  Each  deepwater  wave 
condition  was  represented  as  a  combination  of  many  multi-directional  wave 
components.  The  amount  of  directional  spreading  was  greatest  for  short  period 
cases  and  gradually  narrowed  as  Tp  increased.  Each  component  was  refracted 
over  straight,  parallel  bottom  contours  with  an  east-west  orientation.  Components 
were  recombined  in  13-ft  depth  to  give  a  shallow  water  significant  wave  height 
and  dominant  direction. 

Wave  climate  in  13-ft  depth  seaward  of  the  harbor  is  summarized  in  Figs.  12 
and  13.  The  distribution  of  wave  periods,  Tp  ,  is  not  shown  because  it  is  basically 
unchanged  in  the  transformation  process.  Wave  directions  are  concentrated 
between  the  south  and  south  southwest. 


Significant  Wave  Height  (m) 


Figure  12.  Harbor  entrance  wave  climate,  H, 
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Objectives  and  Approach 

The  numerical  model  study  had  two  objectives: 

a.  Advance  understanding  of  die  existing  harbor  wave  response. 

b.  Evaluate  die  effect  of  proposed  harbor  modifications  on  wave  response. 

The  harbor  wave  response  model  is  presented  in  die  following  section, 
including  a  brief  description  of  die  HARBD  model  and  implementation  of  die 
model  at  Kildaola  Harbor.  The  final  section  of  this  chapter  describes  die  test 
procedures  and  calculations. 

As  part  of  the  test  procedures,  a  suite  of  incident  wave  conditions  must  be 
specified  at  the  seaward  boundary  of  the  area  covered  by  HARBD.  Incident  short 
waves  are  determined  by  consideration  of  transformed  WIS  information  outside 
the  harbor.  Incident  long  waves  are  specified  over  a  broad  range  of  frequencies 
but  only  a  normally-incident  direction  to  identify  possible  harbor  resonant 
responses. 


Model  Description 

Model  Formulation 

The  numerical  wave  model  HARBD  is  a  steady  state  hybrid  element  model 
used  in  die  calculation  of  linear  wave  response  in  harbors  of  varying  size  and 
depth  (Chen  1986,  Chen  and  Houston  1987,  and  Lillycrop  and  Thompson  1996). 
The  model  as  applied  in  this  study  is  described  in  a  recent  report  on  Kahului 
Harbor,  Maui,  Hawaii  (Thompson  et  al.  1996).  An  overview  of  the  model  and  its 
applications  is  given  by  Thompson  and  Hadley  (1995). 

The  principal  output  information  available  from  HARBD  consists  of 
amplification  factor  and  phase  at  each  node  in  die  numerical  grid.  These  are 
defined  as 
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where  =  amplification  factor, 

a,at  =  local  and  incident  wave  amplitudes, 
H,  Hi  =  local  and  incident  wave  heights, 

<t>  =  velocity  potential, 

i/r  =  phase  relative  to  die  incident  wave, 

Im{4> }  =  imaginary  part  of  <p, 

Re{<f>}  =  real  part  of  <f>. 


Amplification  factors  are  easily  interpreted.  Phases  are  helpful  in  viewing  wind 
wave  and  swell  propagation  characteristics  and  in  interpreting  standing  wave 
patterns.  In  long  wave  applications,  phases  prove  useful  for  determining  relative 
phase  differences  within  the  harbor,  interpreting  harbor  oscillation  patterns,  and 
identifying  potentially  troublesome  nodal  areas. 


Spectral  Adaptation 


HARBD  computes  harbor  response  to  specified  wave  period  and  direction 
combinations.  However  the  model  is  often  used  to  approximate  irregular  wind 
wave  and  swell  behavior,  as  in  physical  model  tests  with  irregular  waves  and  all 
field  cases.  More  realistic  numerical  model  simulations  can  be  obtained  by 
linearly  combining  HARBD  results  from  a  range  of  regular  wave  frequencies  and 
directions  in  the  irregular  wave  spectrum.  With  proper  weighting,  regular  wave 
results  represent  a  desired  spectral  distribution  of  energy. 

Spectral  adaptation  of  the  HARBD  model  is  done  as  a  post-processing  step 
using  tiie  standard,  regular  wave  output  from  the  model.  For  a  given  set  of 
incident  wave  directions  representing  the  range  of  possible  approach  directions, 
HARBD  is  run  for  a  number  of  wave  periods  spread  between  the  shortest  period 
consistent  with  grid  resolution  constraints  and  the  longest  swell  period  of  interest. 
Details  of  the  procedure  are  given  by  Thompson  et  al.  (19%). 

The  effective  amplification  factor  at  each  node  is  computed  as 


nt  nd 


EE  W„^c  £np(fk,  o„) 


(2) 


where  (A^)^  =  effective,  or  spectral,  amplification  factor  at  a  node 

^anp(fk>  @n)  =  nodal  amplification  factor  for  HARBD  computational 
frequency  fk  and  direction  6n 

Nt  =  number  of  HARBD  computational  wave  periods 

Nd  =  number  of  HARBD  computational  wave  directions 
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H- 


=  weighting  factor  for  Jt’th  HARBD  computational  frequency 
=  weighting  factor  for  n’th  HARBD  computational  direction 


Finite  Element  Grids 

The  finite  element  numerical  grid  depicting  existing  conditions  at  Kildaola 
Harbor  was  created  using  WES's  finite  element  grid  development  software 
(Turner  and  Baptista  1993)  (Figure  14).  The  grid  covers  die  entire  Kildaola 
Harbor  area  and  extends  somewhat  seaward  into  Waimea  Bay.  The  land 
boundary  was  matched  to  recent  POD  surveys  of  die  harbor.  Grid  element  size  is 
based  on  die  criterion  of  6  elements  per  wavelength  (die  minimum  recommended 
resolution  with  HARBD)  for  a  6-sec  wave  in  4-fit  water  depth.  Depths  for  areas 
of  interest  in  die  existing  and  plan  harbors  generally  exceed  4  ft.  Some  areas  of 
interest  in  die  existing  harbor  are  shallower  than  4  ft  (even  at  a  high  tide  water 
level  of  + 1  ft  MLLW),  but  it  was  impractical  to  make  the  grid  significandy  finer. 
The  grid  was  expected  to  be  adequate  for  die  existing  harbor,  as  well  as  die  plan 
harbors.  For  die  longer  period  waves,  die  grid  gives  a  high  degree  of  resolution. 
Grid  characteristics  are  summarized  in  Table  3. 


Figure  14.  Grid  of  existing  harbor 
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Table  3 

Grid  Sizes 

Harbor 

Plan 

Number  of: 

Length 
of  Typical 
Element 
(ft) 

Elements 

Nodes 

Solid 

Boundary 

Elements 

Semicircle 

Boundary 

Nodes 

Existing 

24227 

12461 

462 

232 

11 

Plan  1 

24524 

12602 

447 

232 

11 

Plan  6 

27345 

14014 

450 

232 

11 

The  radius  of  die  seaward  semicircle  is  854  ft.  This  is  equivalent  to  2. 1  and 
8. 1  wavelengths  for  die  longest  and  shortest  short  wave  periods  considered,  as¬ 
suming  a  representative  water  depth  of  11  ft  (10-ft  depth  below  MLLW  plus  1-ft 
tide).  The  semicircle  size  and  location  were  chosen  to  include  both  breakwaters 
and  the  immediate  nearshore  area.  The  semicircle  extends  sufficientiy  far 
seaward  to  cover  die  most  important  nearshore  bathymetry  while  maintaining  a 
reasonable  number  of  grid  elements. 

Bathymetric  data  were  obtained  from  recent  (summmer  of  1996)  POD  surveys 
of  the  harbor  area  and  extending  seaward  to  die  300  ft  depth  contour.  NOAA 
hydrographic  chart  19386  provided  a  useful  reference  for  bathymetry  outside  die 
survey  area.  Digitized  depths  were  transferred  onto  the  finite  element  grid  using 
the  WES  grid  software  package.  A  contour  plot  of  bathymetry  in  die  existing 
harbor  is  given  in  Figure  15.  Bathymetry  was  modified  for  Plans  1  and  6  to 
include  project  deptiis  in  the  entrance  and  access  channels  (Appendix  A).  Per 
discussions  with  POD,  die  plan  bathymetry  also  includes  deepening  of  the  harbor 
areas  expected  to  be  used  for  mooring  (Fig.  16). 

Reflection  coefficients,  Kr,  are  needed  for  all  solid  boundaries.  For  die  short 
wave  tests,  Kr  values  were  estimated  from  existing  Corps  of  Engineers  guidance, 
photos,  and  field  notes  from  a  May  1996  site  visit  by  WES  and  POD 
representatives,  and  past  experience.  The  solid  boundary  was  divided  into  nine 
zones  and  a  reflection  coefficient  was  estimated  for  each  zone  (Figure  17). 
Reflection  coefficients  range  from  0.05  for  the  shallow,  gendy  sloped  beach  along 
the  southeast  shore  of  the  existing  harbor  to  0.5  for  all  breakwater  areas.  Similar 
reflection  coefficients  were  used  in  die  plan  harbors.  Additional  parameter  values 
used  in  die  numerical  model  are  summarized  in  Table  4. 

Different  parameters  are  used  for  the  long  wave  tests.  The  reflection 
coefficient  was  set  to  1.0  for  all  solid  boundaries,  since  long  waves  generally 
reflect  very  well  from  a  coastal  boundary.  A  reflection  coefficient  of  zero  was 
used  along  die  open  boundary  west  of  die  west  breakwater.  Long  waves  are  more 
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Figure  15.  Bathymetry,  existing  harbor 


■  —  1  — . - 

Table  4 

Parameter  Values  Used  In  HARBD 

Parameter 

Value  H 

BW 

Bottom  friction,  fi 

0.0 

0.032 

Coastline  reflection, 

1.0 

1.0 

10  ft  below  MLLW 

10  ft  below  MLLW 

affected  by  bottom  friction  than  short  waves,  so  a  value  of  P  greater  than  zero  is 
appropriate.  A  value  of  P =0.032  was  selected,  based  on  experience  calibrating 
to  field  data  at  other  sites,  principally  Kahului  Harbor  (Thompson  et  al.  19%). 
This  and  other  long  wave  parameters  are  given  in  Table  4. 
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Figure  17.  Wave  reflection  coefficient  values,  short  waves,  existing  harbor 


Test  Procedures  and 
Calculations 

Incident  Wave  Conditions 

A  range  of  short  and  long  wave 
conditions  incident  to  Kiidaola  Harbor 
was  considered.  A  representative 
range  of  wave  periods  and  directions 
which  could  cause  damaging 
conditions  inside  the  harbor  was 
included,  based  on  WIS  information. 

The  short  wave  periods  and 
approach  directions  considered  are 
given  in  Table  5.  These  conditions 
provide  reasonable  coverage  of  the 
WIS  information  for  the  area, 
summarized  in  Figures  10,  12  and  13. 
The  shortest  period  represents  a  local 


Table  5  1 

Summary  of  Incident  Short  B 
Wave  Conditions  || 

Wave  Period 
(sec) 

Wave  Direction  H 

(deg,  coming  from)  H 

6 

15 

164 

7 

16 

184 

8 

17 

204 

9 

18 

10 

19 

11 

20 

12 

.21 

13 

22 

14 
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sea  condition  and  die  longest  represents  a  very  long  swell  condition.  Directions 
were  chosen  to  include  likely  approach  directions  to  die  harbor  entrance  and  to 
give  adequate  representation  of  the  directional  spectrum  in  post-processing.  Test 
directions  were  straight  into  die  harbor  (184  deg)  and  20-deg  increments  to  either 
side  (coming  from,  relative  to  true  north).  Incident  wave  directions  and  die 
angular  orientation  of  die  seaward  semicircular  model  boundary  are  illustrated  in 
Figure  18. 

For  die  study  of  existing  harbor  conditions  and  comparison  of  alternatives, 
HARBD  was  run  with  the  full  set  of  short  wave  periods  and  directions  in  all 
possible  combinations.  Model  results  were  then  evaluated  for  directional  spectra 

with  Tp  and  6m  values  as  follows:  periods  of  6,  7 . 20  sec  and  directions  of 

•••>  210  deg  azimuth  (coming  from).  These  values  cover  die  range  of 
conditions  in  the  WIS  nearshore  information. 

Incident  long  wave  conditions  considered  are  given  in  Table  6.  A  fine 
resolution  in  wave  frequency  was  used  over  the  full  range  of  possible  resonant 
conditions  to  insure  that  all  important  peaks  were  identified.  A  total  of  451 
periods  was  considered.  Only  one  approach  direction  is  included,  since  past 
studies  have  indicated  that  harbor  response  is  relatively  insensitive  to  incident  long 
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wave  direction.  Tins  direction  represents 
a  wave  directly  approaching  the  harbor 
entrance  from  deep  water. 

One  water  level  was  tested.  The  tide 
range  at  Kikiaola  Harbor  is  relatively 
small,  with  a  mean  range  of  about  1  ft. 
Harbor  wave  response  is  unlikely  to  vary 
much  with  water  level  over  tiiis  tidal 
range.  The  water  level  was  selected  as 
+ 1.0  ft  above  mean  lower  low  water, 
representative  of  a  high  tide  condition. 


Calculation  of  Spectra 

Numerical  model  test  results  for  short 
waves  in  Kikiaola  Harbor  are  all  based  on 
spectral  post-processing  of  the  initial 
HARBD  runs.  Hence,  short  wave 
amplification  factors  are  all  in  the  form  of 
q.  2.  This  approach  requires, 
first,  that  HARBD  be  run  with  the  range 
of  wave  periods  and  directions  to  be 
considered  in  the  spectral  calculations. 

Second,  values  of  peak  wave  period,  Tp,  corresponding  to  the  peak  spectral 
frequency;  wave  approach  direction,  6m;  spectral  peak  enhancement  factor,  y; 
and  directional  spreading  factor,  s,  must  be  specified.  The  Tp  and  6m  values  were 
chosen  to  represent  wind  wave  and  swell  conditions  at  the  harbor,  as  discussed  in 
tiie  section  on  Incident  Wave  Conditions.  Values  for  y  and  s  were  approximated 
as  discussed  by  Thompson  et  al.  (1996). 

Output  Basins 

In  order  to  get  special  coverage  of  areas  where  harbor  traffic  would  most 
likely  be  affected  by  wave  conditions,  between  15  and  25  output  locations  or 
"basins"  were  selected  to  cover  each  harbor  layout.  A  basin  is  a  small  cluster  of 
elements  over  which  the  HARBD  response  is  averaged  to  give  a  more 
representative  output.  Whenever  possible  basins  were  positioned  to  coincide  with 
basins  of  other  plans,  particularly  those  of  the  existing  harbor  (Appendix  A). 

Each  basin  in  tiiis  study  contains  12-15  elements.  HARBD  output  information 
was  saved  at  each  of  these  locations  in  addition  to  the  detailed  output  at  nodes. 


1  Table  6 

Summary  of 

Incident  Long  Wave 
Conditions 

Wave 

Wave  Direction  1 

Period 

(deg,  coming  I 

(sec) 

from)  1 

25.00 

184  1 

25  06 

I 

2513 

| 

i 

I 

500  0 

1 

1  Frequency  Increments  are  H 

0.0001  H2  for  periods  of  25-80  1 
sec  and  0.00006  Hz  for  I 

periods  of  80-500  sec  || 
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4  Harbor  Response  to  Wind 
Waves  and  Swell 


Numerical  model  studies  of  the  harbor  response  to  wind  waves  and  swell  were 
directed  primarily  toward  assessing  the  operational  performance  of  alternative 
harbor  modifications.  Results  are  summarized  in  this  chapter.  Amplification 
factors  are  presented  in  the  following  section.  The  final  section  gives  Hs  values 
exceeded  10  percent  and  1  percent  of  the  time,  a  result  more  directly  applicable  to 
operational  performance.  The  Hs  values  are  derived  from  a  combination  of 
amplification  factors  from  the  numerical  model  and  wave  hindcast  information 
outside  the  harbor.  They  are  compared  to  operational  criteria  for  wind  waves  and 
swell. 

Amplification  Factors 


Amplification  factors,  representing  directionally-spread  short  wave  spectra  in 
the  form  of  (Aanp)e^  in  Eq.  2,  were  calculated  for  a  variety  of  wind  wave  and  swell 
conditions.  Figure  19  illustrates  the  behavior  of  a  common  wave  condition  at 
Kikiaola  Harbor,  a  12-sec  wave  approaching  from  200  deg  azimuth.  Contour 
plots  of  for  the  existing  and  plan  harbors  are  shown.  The  plots  indicate 

reduced  wave  heights  in  the  entrance  channel  near  tiie  main  breakwater  entrance 
in  the  plan  conditions,  as  compared  to  the  shoaled  entrance  in  the  existing  harbor. 
Also,  the  plans  are  more  effective  than  the  existing  harbor  in  affording  protection 
to  wharf  and  mooring  areas  from  this  wave  condition. 

Plots  of  wave  phase  (^in  Eq.  1)  are  included  in  Figure  19.  Since  the  phase 
lines  show  the  alignment  of  wave  crests,  they  give  a  visual  representation  of 
diffraction  and  shoaling  effects  on  wave  direction  and  length  as  the  12-sec  waves 
interact  with  harbor  structures  and  bathymetry.  For  clarity  of  presentation,  the 
phase  plots  were  taken  from  12-sec  regular  wave  results. 

For  a  more  concise  comparison  between  the  existing  harbor  and  alternative 
plans,  average  values  of  were  computed  for  each  basin  across  wave 

periods  ranging  from  6  sec  through  20  sec.  Figure  20  illustrates  results  for  the 
existing  harbor  and  two  plans.  The  changes  progressively  as  incident 

wave  direction  changes.  As  would  be  expected,  amplification  tends  to  be  greater 
for  directions  of  more  direct  approach  to  the  basins.  As  illustrated  in  Figure  19, 
the  wind  wave  and  swell  response  in  the  harbor  is  basically  a  result  of  diffraction 
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Figure  19.  Amplification  factor  and  phase  contours,  12-sec  wave  period,  200-deg  direction 
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Amplification  Factor 


through  the  breakwater  gap.  Boundary  reflection  characteristics  have  only  a 
localized  effect  on  die  waves.  Complete  results  of  (A^J^  averaged  over  wave 
period  are  given  in  Appendix  B. 

An  even  more  concise  description  of  (A^^  at  each  basin  can  be  obtained  by 
considering  wave  climate  as  well.  A  climate-based  amplification  factor  is 
calculated  for  each  basin  as 

*d  V 

climate  ~  £  eff  )i  (3) 


where 

jk  =  indices  denoting  th tf  period  interval  and  indirection  interval,  where  die 
intervals  are  based  on  die  incident  wave  conditions  in  Table  5 
((A =  spectral  amplification  factor  for  die/*  period  and  indirection 
=  number  of  incident  wave  conditions  with  Tp  and  6m  in  the/*  and  k* 
period  and  direction  intervals 
N«<*j  =  total  number  of  incident  wave  conditions 

This  climate-based  amplification  factor  is  given  in  Appendix  B  for  every  basin  and 
harbor  plan. 

Amplification  factors  for  basins  in  die  shallower  harbor  areas  can  be  expected 
to  be  conservative  because  bottom  friction  was  set  to  zero  in  die  HARBD  model 
runs.  This  choice  of  bottom  friction  is  standard  procedure  for  wind  wave  and 
swell  simulations  and  it  has  served  well  in  many  previous  studies.  However,  in 
die  study  of  Kahului  Harbor,  it  appeared  clear  that  die  lack  of  bottom  friction  was 
significantly  affecting  model  results  in  one  very  shallow  basin  (Thompson  et  al. 
1996).  Some  very  shallow  basins  in  Kildaola  Harbor  may  be  similarly  affected. 

A  trial  run  of  the  existing  harbor  with  /?=0.032  indicated  die  extra  complications 
of  using  bottom  friction  (which  makes  amplification  factors  dependent  on  wave 
incident  wave  height  as  well  as  period  and  direction)  were  not  warranted  in  this 
study. 


Evaluation  Against  Operational  Criteria  for  Wind 
Waves  and  Swell 


Standard  operational  criteria  used  by  the  U.S.  Army  Corps  of  Engineers 
(USACE)  for  wind  waves  and  swell  in  shallow  draft  harbors  are: 

-  wave  height  in  berthing  areas  will  not  exceed  1  ft  more  than  10  percent  of  the 
time 

-  wave  height  in  entrance  and  access  channels  and  turning  basins  will  not 
exceed  2  ft  more  than  10  percent  of  the  time 
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Wave  heights  for  assessing  the  USACE  criteria  were  computed  by  combining 
die  time  history  of  wave  hindcast  information  with  harbor  model  results  to  create 
a  time  history  of  wave  heights  at  each  harbor  basin.  For  each  WIS  hindcast  time, 
die  corresponding  wave  height  at  a  harbor  basin  is 


(fy  harbor  ~  (4imp)qf  X  (fywiS  (4) 


where 

(HJharbor  —  significant  wave  height  at  a  harbor  basin 

—  spectral  amplification  factor  interpolated  between  values  for 
periods  and  directions  in  Table  5  to  represent  Tp  and  6m  at  the  seaward 
HARBD  boundary 

(HJms  =  significant  wave  height  at  the  seaward  HARBD  boundary 


The  12-month  time  history  of  (HJharbor  at  each  basin  was  sorted  into  descending 
order  and  die  value  of  H,  which  was  exceeded  10  percent  of  the  time  was 
identified.  The  Hs  value  exceeded  1  percent  of  die  time  was  also  identified.  The 
H}  with  1  percent  exceedance  relates  to  a  more  demanding  operational  condition, 
which  may  also  be  helpful  in  assessing  performance  of  the  harbor  plans. 

Significant  wave  heights  exceeded  10  percent  of  the  time  are  less  than  1  ft  at 
die  wharf  in  die  existing  and  plan  harbors  (Figure  21).  The  USACE  criteria  for 
acceptable  berthing  areas  and  channels  are  shown  in  the  figure  as  solid  lines. 

Wave  conditions  in  die  inner  channel  satisfy  die  USACE  criterion  in  the  existing 
and  plan  harbors.  However,  wave  conditions  in  die  outer  portion  of  the  existing 
channel  exceed  the  criterion.  Wave  conditions  in  the  outermost  portions  of  the 
channel  in  Plans  1  and  6  slighdy  exceed  the  USACE  criterion.  In  conjunction 
with  die  increased  width  of  die  outermost  portion  of  die  plan  channels,  the  small 
exceedance  of  the  USACE  criterion  is  unlikely  to  interfere  with  safe  navigation. 

The  Hs  values  exceeded  1  percent  of  die  time  are  considerably  higher  than 
those  exceeded  10  percent  of  the  time,  but  show  similar  relative  trends  (Appendix 
B).  The  existing  wharf  area  still  falls  below  the  1-ft  wave  height  threshold. 

Values  of  Hs  with  10-  and  1-percent  exceedence  are  tabulated  in  Appendix  B. 
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Figure  21 .  Comparison  of  Ht  exceeded  1 0  percent  of  the  time 
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5  Harbor  Oscillations 


To  evaluate  harbor  resonance  characteristics,  the  HARBD  numerical  model 
was  run  for  the  existing  harbor  and  Plan  6.  Incident  long  wave  periods  ranged 
from  25  sec  to  495  sec  in  very  fine  increments,  as  discussed  in  Chapter  3.  These 
evaluations  were  included  because  oscillations  are  an  important  part  of 
interpreting  die  existing  harbor  wave  response,  and  modifications  to  the  harbor 
can  potentially  lead  to  operational  problems  due  to  harbor  oscillations. 
Amplification  factor  results  are  presented  in  the  following  section.  The 
operational  significance  of  the  results  is  discussed  in  the  final  section. 

Amplification  Factors 


Amplification  factors  for  the  long  waves 
involved  in  harbor  oscillation  behave 
differently  than  those  for  wind  waves  and 
swell.  Long  waves,  because  of  their  length 
relative  to  harbor  dimensions  and  their 
reflectivity  from  harbor  boundaries,  form 
standing  wave  patterns  in  the  harbor. 

Standing  wave  behavior  in  a  simple  closed 
basin  of  uniform  depth  is  illustrated  in  Figure 
22.  In  tiie  fundamental  mode  of  oscillation, 
antinodes  occur  at  both  basin  walls  and  a  node 
midway  between  walls.  Second  and  third 
modes  of  oscillation  are  also  illustrated. 
Antinodes  always  occur  at  the  walls. 
Additional  antinodes  and  nodes  occur  at 
regular  intervals  between  walls,  with  the 
number  of  antinodes  and  nodes  dependent  on 
the  mode  of  oscillation. 


Figure  22.  Harbor  oscillation 
definitions 


The  water  surface  in  a  standing  wave  has 
its  greatest  vertical  motion  at  antinodes. 
There  is  no  vertical  movement  at  an  ideal 
node,  but  horizontal  velocities  reach  a 
maximum  there.  In  terms  of  amplification 
factors,  this  behavior  gives  large  values  of 
at  antinodes  and  small  values  around 
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nodes  (die  notation  A^  ,  is  used  to  denote  long  wave  amplification  factors). 
Contrary  to  wind  waves  and  swell,  small  values  of  A^,  are  not  necessarily 
indicative  of  a  tranquil  harbor  area. 

Phases  in  a  standing  wave  also  behave  differently  than  for  typical  wind  waves 
and  swell.  For  example,  die  water  surface  in  die  fundamental  mode  of  oscillation 
in  Figure  22  simultaneously  reaches  a  maximum  at  every  point  to  die  left  of  the 
node.  These  points  are  all  in  phase.  At  the  same  time,  every  point  to  die  right  of 
die  node  reaches  a  minimum  value.  These  points  are  also  in  phase  with  each 
other  but  exactiy  out  of  phase  with  the  points  to  die  left  of  die  node.  Thus  phases 
in  a  simple  standing  wave  are  constant  between  an  antinode  and  node.  They 
quickly  change  by  180  deg  (or  ;r  radians)  across  the  node  and  remain  constant  up 
to  die  next  node  or  boundary. 

Because  of  the  nature  of  long  wave  behavior,  harbor  oscillation  studies  must 
include  consideration  of  likely  antinode  locations,  particularly  harbor  corners. 
Several  output  basins  were  added  to  those  used  for  wind  waves  and  swell,  as 
shown  in  Appendix  A. 

Amplification  factors  for  basins  in  the  existing  harbor  and  Plan  6  are  shown  as 
a  function  of  wave  frequency  in  Appendix  C.  Some  frequencies  produce  a  strong 
resonant  amplification,  with  peak  amplification  factors  between  about  2  and  7. 
Some  of  die  same  resonant  frequencies  appear  at  several  basins  though  die 
strength  of  amplification  can  vary  considerably  between  basins.  A  large  peak  at 
very  low  frequency  (0.00214  Hz  or  470-sec  period  in  existing  harbor;  0.00274  Hz 
or  365-sec  period  in  Plan  6)  shows  at  every  basin  and  plan.  This  peak  represents 
die  Helmholtz  (or  grave)  mode  of  oscillation,  in  which  die  entire  harbor  rises  and 
falls  in  unison.  Phase  is  constant  over  the  whole  harbor.  Computed  values  of 
Atm^.i  over  the  range  of  frequencies  up  to  and  including  die  Helmholtz  mode  were 
divided  by  two  because  these  oscillations  affect  die  entire  numerical  model  domain 
and  would  otherwise  give  a  distorted  view  of  die  harbor  effect. 

Amplification  factor  and  phase  contour  plots  for  six  of  the  dominant  resonant 
peaks  in  die  existing  harbor  (excluding  Helmholtz  resonance)  show  oscillation 
patterns.  In  die  amplification  factor  plots,  areas  of  high  amplification  are  evident 
as  darker  shades  of  gray  (Figure  23).  Corresponding  phase  contours  are  shown  in 
Figure  24.  Areas  in  which  phase  contours  are  tightly  bunched  indicate  nodal 
areas.  As  would  be  expected  for  standing  waves,  nodal  lines  in  Figure  24 
coincide  with  low  amplification  factors  in  Figure  23.  The  phase  plots  also 
indicate  areas  of  the  harbor  which  rise  and  fall  together  during  die  resonant 
condition  (same  gray  shade).  Thus  die  oscillation  patterns  can  be  interpreted. 

The  150.60-sec  resonant  period  shown  in  Figure  23  represents  a  relatively 
simple  rocking  oscillation  between  the  outer  and  inner  harbor  areas.  The  nodal 
line  lies  just  inside  of  the  entrance  to  die  inner  harbor.  In  die  61 .35-sec  period 
resonance,  opposite  comers  of  the  harbor  complex  act  in  phase.  The  main 
entrance  and  northeast  comer  rise  and  fall  together,  as  do  die  north  part  of  the 
outer  harbor  and  the  southeast  comer  of  the  inner  harbor.  A  nodal  line  intersects 
die  wharf.  Shorter  period  resonances  occur  across  shorter  harbor  dimensions  or 
as  higher  order  modes  along  die  longer  dimensions.  The  40. 16-sec  resonance 
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Figure  23.  Resonant  long  wave  amplification  factor  contours,  existing  harbor 
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Figure  24.  Resonant  long  wave  phase  contours,  existing  harbor 
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between  die  west  breakwater  and  the  inner  east  breakwater  stub  is  of  special  note 
because  it  generated  die  highest  amplification  factors. 

Long  wave  amplification  factors  shown  here  may  be  overestimated  for 
resonant  peaks  at  periods  less  than  about  100  sec.  Wave  reflection  coefficient  at 
all  solid  boundaries  was  set  to  1.0  for  all  long  wave  runs,  but  die  recent  Kahului 
Harbor  study  and  physical  reasoning  indicate  that  small  reductions  in  reflection 
coefficient  at  the  shorter  long  wave  periods  may  be  appropriate. 

Amplification  factor  and  phase  contour  plots  for  die  main  resonant  frequencies 
in  Plan  6  are  given  in  Figures  25  and  26.  The  longest  period  resonance,  with  a 
period  of  113.64  sec,  is  a  simple  rocking  motion  between  die  inner  and  outer 
harbors.  The  period  of  this  resonance  is  shorter  than  the  similar  resonance  in  the 
existing  harbor,  mainly  due  to  the  greater  water  depths  in  Plan  6.  The  73.53-sec 
resonance  is  a  rocking  motion  between  die  east  and  west  areas  of  the  inner 
harbor.  The  north  area  of  the  outer  harbor  also  participates  in  phase  with  die  east 
inner  harbor.  The  57.47-sec  motion  is  a  simple  rocking  between  the  north  and 
south  areas  of  die  harbor  complex.  As  with  die  existing  harbor,  a  strong 
resonance  occurs  between  die  west  breakwater  and  die  inner  east  breakwater  stub. 
The  period  of  this  resonance,  39.68  sec,  is  nearly  die  same  as  in  die  existing 
harbor. 


Evaluation  Against  Operational  Criteria  for  Long 
Waves 


Procedures  for  evaluating  die  operational  acceptability  of  different  harbor 
plans  subjected  to  long  waves  are  reviewed  by  Thompson  et  al.  (1996)  in  relation 
to  die  deep  draft  Kahului  Harbor.  Long  wave  heights  (which  are  unavailable  at 
Kikiaola  Harbor  since  no  wave  gage  data  were  collected)  are  a  key  factor  in  most 
of  the  procedures,  but  some  operational  evaluation  can  be  based  on  amplification 
factors. 

An  operational  guideline  is  based  on  the  value  of  Arrr.f  for  the  higher  resonant 
peaks.  Experience  with  Los  Angeles  and  Long  Beach  harbors  has  indicated  that  if 
^anp,i  is  greater  than  about  5,  some  operational  difficulties  may  be  encountered.  If 
A^i  is  greater  than  10,  major  operational  problems  can  be  expected.1  This 
guideline  may  be  applied  to  the  plots  in  Appendix  C.  The  only  output  basins 
signaling  potential  operational  problems  are  located  in  die  comer  just  west  of  the 
inner  east  stub  breakwater  (Basin  21  in  the  existing  harbor;  Basin  24  in  Plan  6). 
Values  of  Alinp  l  at  this  location  are  approximately  7  in  die  existing  harbor  and  5  in 
Plan  6.  Since  this  area  is  not  planned  for  operational  use,  the  existing  and  Plan  6 
harbors  do  not  appear  to  suffer  from  detrimental  harbor  oscillations. 


1  Personal  Communication,  William  C.  Seabergh,  Research  Hydraulic  Engineer,  U.S.  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  MS. 
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Figure  26.  Resonant  long  wave  phase  contours.  Plan  6 
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Plan  6  shows  a  small  improvement  in  harbor  oscillation  performance  over  the 
existing  harbor.  The  differences  are  most  notable  as  reduced  amplification  factor 
peaks  in  harbor  and  outer  basin  comers  (Appendix  C). 
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6  Conclusions  and 
Recommendations 


Numerical  model  studies  of  the  wave  response  of  Kikiaola  Harbor  have 
produced  valuable  information  about  the  existing  harbor  and  possible 
modifications.  The  numerical  model  was  used  to  simulate  the  behavior  of  the 
existing  harbor  and  two  alternative  modifications  to  the  harbor.  Model  results  are 
compared  with  criteria  for  operational  acceptability  and  with  experience  in  the 
existing  harbor  to  the  extent  possible.  The  effectiveness  of  proposed  alternatives 
for  wind  wave  and  swell  protection  does  not  necessarily  correlate  with  protection 
from  oscillations.  These  two  aspects  of  harbor  operability  are  both  considered  in 
judging  success  of  the  alternative  plans. 

Performance  of  the  alternative  plans  can  be  summarized  by  their  success 
relative  to  simple,  meaningful  criteria.  For  wind  waves  and  swell,  success  was 
defined  as  meeting  or  bettering  the  USACE  criteria  for  harbor  and  entrance 
channel  tranquillity  (Chapter  4).  A  major  reduction  in  the  occurrence  of  wave 
breaking  in  the  entrance  channel  is  also  desired.  Plan  performance  relative  to 
each  of  these  wind  wave  and  swell  concerns  is  as  follows: 

a.  USACE  criterion  in  berthing  areas.  Plan  1,  Plan  6,  and  the  existing  harbor 
all  satisfied  this  criterion.  Berthing  areas  are  well  protected  from  wind  waves 
and  swell  in  all  plans. 

b.  USACE  criterion  in  entrance  and  access  channels.  Plan  1  satisfied  this 
criterion  everywhere  except  at  the  seaward  end  of  the  entrance  channel,  where 
the  10  percent  exceedance  Hs  is  2.2  ft  (compared  to  2.0  ft  specified  in  the 
criterion).  Plan  6  also  satisfied  the  criterion  everywhere  except  in  the  seaward 
part  of  the  entrance  channel.  The  existing,  shoaled  outer  entrance  channel 
significantly  failed  the  USACE  criterion.  The  inner  portion  of  the  existing 
channel,  beginning  at  around  the  90-deg  turn,  satisfied  the  criterion. 

For  both  Plan  1  and  Plan  6,  the  outer  entrance  channel  flares  out  to  give  extra 
width  for  vessel  maneuvering.  A  small  exceedance  of  the  USACE  criterion  in 
the  wider,  more  open  channel  may  not  be  cause  for  concern. 

Wave  conditions  in  the  landward  portion  of  the  Plan  6  entrance  channel  are  a 
potential  navigation  concern.  From  where  the  channel  passes  just  west  of  the 
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main  breakwater  up  to  the  90-deg  turn  into  the  lee  of  the  breakwater,  wave 
conditions  approach  (but  do  not  exceed)  die  USACE  criterion.  The  ability  of 
vessels  to  navigate  safely  past  die  breakwater  head  and  turn  into  protected 
waters  should  be  given  due  consideration.  Experience  with  the  existing 
harbor,  which  also  requires  a  90-deg  turn  in  an  area  where  wave  heights  can 
be  fairly  high  (Figure  21),  should  be  helpful  in  assessing  die  navigation 
concerns. 

c.  Breaking  waves.  Breaking  waves  are  a  potential  concern  in  die  exposed 
entrance  channel.  The  HARBD  numerical  model  does  not  identify  breaking 
wave  conditions,  but  they  can  be  inferred  from  wave  heights  and  water  depths. 
In  die  existing  channel,  die  shoaled  depths  are  less  than  6  ft  in  some  areas. 

The  10  percent  exceedance  H,  outside  the  breakwater  entrance  was  estimated 
as  3.6  ft  (Figure  21).  When  H,  exceeds  about  0.5-0.6  times  the  depth, 
signficant  wave  breaking  can  be  expected.  Thus  die  existing  entrance  channel 
would  be  expected  to  experience  breaking  wave  conditions  more  than  10 
percent  of  die  time. 

In  Plan  1  and  Plan  6,  die  channel  depth  is  12  ft  and  die  10  percent  exceedance 
H,  is  around  2  ft  or  less.  Thus  breaking  wave  conditions  would  not  be 
expected  in  these  plans. 

There  are  several  key  limitations  on  the  above  conclusions.  First,  die  coast 
around  Kikiaola  Harbor  is  an  area  of  active  sediment  movement.  The  existing 
channel  and  harbor  have  experienced  shoaling  problems.  The  present  study  was 
based  on  actual  depths  of  the  existing  harbor  but  project  depths  of  Plan  1  and  Plan 
6.  If  either  plan  were  constructed  and  significant  shoaling  occurred,  die  wind 
wave  and  swell  response  would  change  from  die  estimates  of  this  study.  Also, 
wave  breaking  could  again  become  a  problem.  Sediment  transport  and  channel 
shoaling  were  outside  die  scope  of  this  study,  but  they  are  important 
considerations  in  selecting  a  final  plan  for  the  harbor. 

Another  limitation  is  the  wave  climate  information.  It  was  derived  from  one 
year  of  hindcasts  in  a  rather  difficult  area.  Conclusions  based  on  wave  climate 
would  have  had  a  higher  confidence  level  if  high-quality  local  measurements  were 
available.  Limitations  associated  with  die  numerical  harbor  model  were  presented 
in  Chapter  1.  The  absence  of  wave  breaking  and  currents  are  the  most  relevant 
limitations  in  the  present  study. 

Harbor  response  to  long  waves  (harbor  oscillations)  was  analyzed  for  die 
existing  harbor  and  Plan  6.  Plan  1  was  not  included  in  this  part  of  die  study. 
Based  on  experience  in  other  harbors  correlating  long  wave  amplification  factors 
with  operational  problems,  die  existing  harbor  and  Plan  6  both  appear  to  be  free 
of  oscillation  problems.  The  only  harbor  area  of  potential  concern  is  between  die 
west  breakwater  and  the  inner  east  breakwater  stub,  which  is  not  planned  as  a 
mooring  facility.  Oscillations  in  Plan  6  are  less  than  in  the  existing  harbor. 

A  physical  model  study  to  optimize  die  preferred  plan  may  be  cost  effective  as 
a  final  phase  of  the  studies.  A  focused  physical  model  study  is  likely  to  result  in 
cost  savings  beyond  die  study  cost  in  fine-tuning  design  details  such  as  breakwater 
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length,  crest  height,  and  slope.  Also,  die  physical  model  can  reproduce  some 
processes  not  represented  in  the  numerical  model,  such  as  wave  breaking. 
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Appendix  A 

Bathymetry  and  Output  Basin 
Locations 


Figure  A1.  Bathymetry,  Plan  1 
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Appendix  B 
Wind  Wave  and  Swell 
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Model 
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Harbor  Oscillation  Summaries 
from  Numerical  Model 


Appendix  C 


ci 


Amplification  Factor  3  Amplificalion  Factor 


'3.  Long  wave  response,  harbor  center,  existing  harbor 


Figure  C4.  Long  wave  response,  harbor  center.  Plan  6 
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Appendix  D 
Notation 


a  Wave  amplitude,  m  (ft) 

Oi  Incident  wave  amplitude,  m  (ft) 

Wave  amplification  fector 

(Aanp)tff  Effective,  or  spectral,  wave  amplification  factor 

A^,  Wave  amplification  fector,  long  waves  (harbor  oscillations) 

d  Water  depth,  m  (ft) 

dp  Water  depth  for  away  from  harbor,  m  (ft) 

f  Wave  frequency,  Hz 

H  Wave  height,  m  (ft) 

Hj  Incident  wave  height,  m  (ft) 

Ht  Significant  wave  height,  m  (ft) 

Kr  Reflection  coefficient  of  a  solid  boundary 

Kr  ceas,  Reflection  coefficient  of  coastline  fer  away  from  harbor 

Nd  Number  of  HARBD  computational  wave  directions  for  spectral  approximation 

Nt  Number  of  HARBD  computational  wave  periods  for  spectral  approximation 

s  Directional  spreading  fector 

Tm  Mean  spectral  wave  period,  sec 

Tp  Peak  spectral  wave  period,  sec 

M>k  Weighting  fector  for  k’th  HARBD  computational  frequency 
P  Dimensionless  bottom  friction  coefficient 
Y  Spectral  peak  enhancement  fector 

6m  Mean  wave  direction 

<p  Velocity  potential 

ifr  Wave  phase 
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